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We analyze characteristics of drifter trajectories from the Adriatic Sea with 
recently introduced nonlinear dynamics techniques. We discuss how in quasi- 
enclosed basins, relative dispersion as function of time, a standard analysis 
tool in this context, may give a distorted picture of the dynamics. We further 
show that useful information may be obtained by using two related non- 
asymptotic indicators, the Finite-Scale Lyapunov Exponent (FSLE) and the 
Lagrangian Structure Function (LSF), which both describe intrinsic physical 
properties at a given scale. We introduce a simple chaotic model for drifter 
motion in this system, and show by comparison with the model that La- 
grangian dispersion is mainly driven by advection at sub-basin scales until 
saturation sets in. 
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I. INTRODUCTION 

Understanding the mechanisms of transport and mixing processes is an important and 
challenging task, which has wide relevance from a theoretical point of view, e.g. for the study 
of diffusion and chaos in geophysical systems in general, or for validating simulation results 
from a general circulation model. It is also a necessary tool in the analysis of problems of 
general interest and social impact, such as the dispersion of nutrients or pollutants in sea 
water with consequent effects on marine life and on the environment (Adler et al., 1996). 

Recently, a number of oceanographic programs have been devoted to the study of the 
surface circulation of the Adriatic Sea by the observation of Lagrangian drifters, within the 
larger framework of drifter-related research in the whole Mediterranean Sea (Poulain, 1999). 

The Adriatic Sea is a quasi-enclosed basin, about 800 long by 200 km wide, connected to 
the rest of the Mediterranean Sea through the Otranto Strait. From a topographic point of 
view, three major regions can be considered: the northern part is the shallowest, about 100 
m maximum depth, and extends down to the latitude of Ancona; the central part, which 
deepens down to about 260 m in the Jabuka Pit, and the southern part which extends from 
the Gargano promontory to the Otranto Strait. The southern part is the deepest, reaching 
about 1200 m in the South Adriatic Pit. Reviews on the oceanography of the Adriatic Sea 
can be found in Artegiani et al. (1997), Orlic et al. (1992), Poulain (1999) and Zore (1956). 

Lagrangian data offer the opportunity to employ techniques of analysis, well established 
in the theory of chaotic dynamical systems, to study the behavior of actual trajectories and 
compare those with a kinematic model. 

Let us assume that the Lagrangian drifters are passively advected in a two-dimensional 
flow, e.g. as would be the case in a frictionless barotropic approximation (Ottino, 1989; 
Crisanti et al., 1991): 

— = u{x,y,t) and — = v{x,y,t) , (1) 

where {x{t),y{t)) is the position of a fluid particle at time t in terms of longitude and latitude 
and u and v are the zonal and meridional velocity fields, respectively. 



For the Eulerian description of a geophysical system, one should in principle use numer- 
ical solutions of the Navier-Stokes equations (or other suitable equations, e.g. the quasi- 
geostrophic model) to obtain the velocity fields. In practice, direct numerical simulation 
of these equations on oceanographic length scales is of course not possible, and one has to 
invoke approximations, i.e. turbulence modeling. This motivates to use instead a simpli- 
fied kinematic approach, by adopting a given Eulerian velocity field. The criteria for the 
construction of such a field follow from phenomenological arguments and/or experimental 
observation, and have recently been reviewed in (Yang, 1996; Samelson, 1996). 

Let us consider the relationship between Eulerian and Lagrangian properties of a system. 
A wide literature on this topic (e.g. Ottino, 1989; Crisanti et al., 1991) allows us state that, 
in general, motion in Eulerian and Lagrangian variables can be rather different. It is not 
rare to have regular Eulerian behavior, e.g. a time-periodic velocity field, co-existing with 
Lagrangian chaos or vice- versa. 

In quasi-enclosed basins like the Adriatic Sea, a characterization of the mechanisms of the 
mixing is highly non-trivial. We first observe, see below for detailed discussion, that the use 
of the standard diffusion coefficients can have rather limitated applicability (see Artale et al., 
1997). Already classical studies on Lagrangian particles in ocean models contain remarks on 
the intrinsic difficulties in using one-particle diffusion statistics (Taylor, 1921). In situations 
where the advective time is not much longer than the typical decorrelation time scale of 
the Lagrangian velocity, the diffusivity parameter related to small-scale turbulent motion 
cannot converge to its asymptotic value (Figueroa and Olson, 1994). One the other hand, 
a generalization of the standard Lyapunov exponent, the Finite-Scale Lyapunov Exponent 
(FSLE), originally introduced for the predictability problem (Aurell et al., 1996:1997), has 
been shown to be a suitable tool to describe non- asymptotic properties of transport. This 
finite-scale approach to Lagrangian transport measures effective rates of particle dispersion 
without assumptions about small-scale turbulent processes. For an alternative method see 
Buffoni et al. (1997); for a recent review and systematic discussion of non- asymptotic 
properties of transport and mixing in realistic cases, see Boffetta et al. (2000). 



In this paper we report data analysis of surface drifter motion in the Adriatic Sea using 
relative dispersion, FSLE and Lagrangian Structure Function (LSF), a quantity related to 
the FSLE. We also introduce a chaotic model for the Lagrangian dynamics, and use the 
FSLE and LSF characteristics to compare model and data. We show that it can be very 
difficult to get an estimate of the diffusion coefficient in a quasi-enclosed basin, and/or to 
look for deviations from the standard diffusion law. In fact, the time a cluster of particles 
takes to spread uniformly and reach the boundaries is not much longer than the largest 
characteristic time of the system. In contrast, the FSLE and the LSF do characterize the 
transport properties of Lagrangian trajectories at a fixed spatial scale. We will finally show 
that a simple kinematic model reproduces the data. 

In section 2 we describe the data set we have used, review relevant concepts and analysis 
techniques for Lagrangian transport and chaos. In section 3 we introduce a kinetic model 
of the Lagrangian dynamics, and in section 4 we compare the data and the model. Section 
5 contains a summary and a discussion of the results. 

II. DATA SET AND ANALYSIS TECHNIQUES 

A. Data set 

In a large drifter research program in Mediterranean Sea, started in the late 80 's and 
continued into the 90's, Lagrangian data from surface drifters deployed in the Adriatic sea 
have been recorded from December 1994 to March 1996. These drifters are similar to the 
CODE (COastal Dynamics Experiment) system (Davis, 1985) and they are designed to 
be sufficiently wind-resistant so to effectively give a description of the circulation at their 
actual depth (1 meter). The drifters were tracked by the Argos Data Location and Collection 
System (DCLS) carried by the NOAA polar-orbiting satellites. It is assumed that after data 
processing drifter positions are accurate to within 200-300 m, and velocities to within 2-3 
cm/s. 



For a description of the experimental program, see Poulain (1999). Technical details 
about the treatment of raw data can also be found in Hansen and Poulain (1996), Poulain 
et al. (1996) and Poulain and Zanasca (1998). 

The data have been stored in separate files, one for each drifter. In the format used by us, 
each file contains: number of records (i.e. number of points of the trajectory); time in days; 
position of the drifter in longitude and latitude; velocity of the drifter along the zonal and 
meridional directions; and temperature in centigrade degrees. The sampling time is 6 hours. 
We can identify five main deployments on which we will concentrate our attention. Selecting 
the tracks by the time of the first record, it is easy to verify that these five subsets consist 
of drifters deployed in the same area in the Otranto strait, near 19 degrees longitude east 
and 40 degrees latitude north. The experimental strategy of simultaneously releasing the 
drifters within a distance of some kilometers, allows us to study dispersion quantitatively. 

From a qualitative point of view, what we observe from the plot of all the trajectories 
(Fig. 1) is the shape of two (cyclonic) basin-wide gyres, located in the middle and southern 
regions respectively, and, an anti-clockwise boundary current which moves the drifters north- 
westward along the east coast and south-eastward down the west coast. The latter is a 
permanent feature of the Adriatic sea (Poulain, 1999). On the other hand, it is known that, 
within a year, the pattern of basin-wide gyres may change between one, two and three gyres 
over a time-scale of months. The southern gyre is the most steady of the three. The data 
also suggest the presence of small scale structures, even though these are likely much more 
variable in time. The time-scale of the typical recirculation period around a basin-wide gyre 
is about one month and the time needed to travel along the coasts and complete one lap of 
the full basin is of the order of a few months. 

B. Analysis techniques 

We recall here some basic concepts about dynamical systems, diffusion and chaos, and 
the quantities that we shall use to characterize the properties of Lagrangian trajectories. 
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If we have Nc clusters of initially close particles, each cluster containing Uk elements, 
relative dispersion can be characterized by the diffusion coefficient 

A = lim,^oo^^^(t) (2) 

with 

sm = ^T.^ E(^f'Ht)- < x,{t) >^'^f (3) 

where 

<a;.(t)>('=)=-E4'"^W (4) 

Xj^ '"' is the i-th spatial coordinate of the j-th particle in the k-th cluster; S"^ = Y.i Sf is the 
mean square displacement of the particles relatively to their time evolving mean position. 
If 5{t) is the distance between two trajectories x*^^^ and x*^^) in a cluster at time t, relative 
dispersion is defined as 

< 5\t) > = < ||x«(t) - x(2)(t)||2 > (5) 

where the average is over all pairs of trajectories in the cluster. In a standard diffusive 
regime, x*^-'^^(t) and x*^^)(t) become independent variables and, for t -^ oo, we have < 5^{t) >= 
2-S'^{t). In the following, we shall consider the cluster mean square radius S'^{t) as a measure 
of relative dispersion. Absolute dispersion, which is defined as the mean square displacement 
from an initial position will not be taken into account in our analysis. 

If, in the asymptotic limit, Sf{t) ~ t^" with a = 1/2, we have the linear law of standard 
diffusion for the mean square displacement, and the Di^s are finite; if a 7^ 1/2 we have 
so-called anomalous diffusion (Bouchaud and Georges, 1990). 

The difficulty that often arises when measuring the exponent a is that, because of the 
finite size of the domain, dispersion cannot reach its true asymptotic behavior. In other 
words, diffusion may not be observable over sufficiently large scales, i.e. much larger than 
the largest Eulerian length scale, and, therefore, we cannot have a robust estimate of the 
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exponent of the asymptotic power law. Moreover, the relevance of asymptotic quantities, 
like the diffusion coefficients, is questionable in the study of realistic cases concerning the 
transport problem in finite-size systems (Artale et al., f997). 

The diffusion coefficients characterize long-time (large-scale) dispersion properties. In 
contrast, at short times (small scales) the relative dispersion is related to the chaotic behavior 
of the Lagrangian trajectories. 

A quantitative measure of instability for the time evolution of a dynamical system (Licht- 
enberg and Lieberman, f992) is commonly given by the Maximum Lyapunov Exponent 
(MLE) A, which gives the rate of exponential separation of two nearby trajectories 

A = hmi^oo hm5(o)_>o - In ^r^ (6) 

where S{t) = ||x*^^)(t) — x'^^^(t)|| is the distance between two trajectories at time t. When A > 
the system is said to be chaotic. There exists a well established algorithm to numerically 
compute the MLE introduced by Benettin et al. (1980). 

A characteristic time, Tx, associated to the MLE is the predictability time, defined as 
the minimum time after which the error on the state of the system becomes larger than a 
tolerance value A, if the initial uncertainty is 6 (Lichtenberg and Lieberman, 1992): 

T. 4 taf (7) 

Let us recall that A is a mathematically well-defined quantity which measures the growth 
of infinitesimal errors. In physical terms, at any time, 6 has to be much lesser than the 
characteristic size of the smallest relevant length of the velocity field. For example, in 3-D 
fully developed turbulence 6 has to be much smaller than the Kolmogorov length. 

When the uncertainty grows up to non-infinitesimal sizes, i.e. macroscopic scales, the 
perturbation 6 is governed by the nonlinear terms and that renders its growth rate a scale- 
dependent index (Aurell et al., 1996, 1997; Artale et al., 1997). It is useful to introduce 
the Finite Scale Lyapunov Exponent (FSLE), X{S). Assuming r > 1 is a fixed amplification 
ratio and < Tr{6) > the mean time that 6 takes to grow up to r ■ 6, we have: 
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The average < ■ > is performed over all the trajectory pairs in a cluster. We note the 
following properties of the FSLE: 

a) in the limit of infinitesimal separation between trajectories, 5-^0, the FLSE tends to 
the maximum Lyapunov exponent (MLE); 

b) in case of standard diffusion, < 5{tY >~ t, we find that \{5) ~ 5^"^ and the propor- 
tionality constant is of the order of the diffusion coefficient; 

c) any slope > —2 for \{5) vs 5 indicates super-diffusive behavior, i.e. non-neglectable 
correlations persist at long times and advection is still relevant; 

d) in particular, when \{5) = constant over a range of scales, we have exponential sep- 
aration between trajectories at constant rate, within that range of scales (chaotic 
advection) . 

Another interesting quantity related to the FSLE is the Lagrangian Structure Function 
(LSF) 1^(6), defined as 

where the value of the velocity difference is taken at the times for which the distance between 
the trajectories enters the scale 6 and the average is performed over a large number of 
realizations. The LSF, z/((5), is a measure of the velocity at which two trajectories depart 
from each other, as a function of scale. By dimensional arguments, we expect that the LSF 
is proportional to the scale of the separation and to the FSLE: 

u{6) ~ 6 X{6) (10) 

so that we should find similar behavior for X{6) and i'{6)/6, if independently measured. 

In order to study the transport properties of the drifter trajectories, we have focused our 
interest on the measurement of Sf{t), X{6) and z^(5). 

9 



With regards to the practical definitions of the FSLE and the LSF, we have chosen a 
range of scales 6 = {6o,6i, ...,6n) separated by a factor r > 1 such that (5j+i = r ■ 6i for 
i = 0,n — 1. The ratio r is often referred to as the "doubling" factor even though it is not 
necessarily equal to 2, e.g. in our case we fixed it at \/2- The r value has naturally an 
inferior bound because of the temporal finite resolution of the trajectories (i.e. it cannot be 
arbitrarily close to 1) and it must be not much larger than 1, if we want to resolve scale 
separation in the system. 

The smallest threshold, 6q, is placed just above the initial mean separation between two 
drifters, ~ 10 km, and the largest one, (5„, is naturally selected by the finite size of the 
domain, ~ 500 km. 

Following the same procedure, it is straightforward to compute the LSF as the mean 
velocity difference between two trajectories at the moment in which the separation reaches 
a scale S: 



jy{S) = < ^{u^-u^Y + iv^-v^Y >s (11) 

where the average is performed over the number of all the pairs within a set of particles, at 
the time in which ||x' — x|| = 6. 

In section 4 below we shall show the results of our data analysis and compare them with 
the simulations from our chaotic model for the Lagrangian dynamics of the Adriatic drifters. 

III. THE CHAOTIC MODEL 

In phenomenological kinetic modeling of geophysical fiows, two possible approaches can 
be considered: stochastic and chaotic. Both procedures generally involve a mean velocity 
field, which gives the motion over large scales, and a perturbation, which describes the 
action of the small scales. The model is stochastic or chaotic if the perturbation is a random 
process or a deterministic time- dependent function, respectively. Examples on kinematic 
mechanisms proposed to model the mixing process can be found in Bower (1991), Samelson 
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(1992), Bower and Lozier (1994), Cencini et al. (1999). 

The choice of one or the other depends on what one is interested in, and what experi- 
mental information is available. In our case, we have opted for a deterministic model since 
there are indications that, at the sea surface, the instabilities of the Eulerian structures are 
mostly due to air-sea interactions, which are nearly periodic perturbations. 

We want to consider a simple model, so let us assume as main features of the surface 
circulation the following elements: an anti-clockwise coastal current; two large cyclonic 
gyres; and some natural irregularities in the Lagrangian motion induced by the small scale 
structures. 

Let us notice that the actual drifters may leave the Adriatic sea through the Otranto 
Strait, but we model our domain with a closed basin, in order to study the effects of the 
finite scales on the transport, and treat it like a,2 — D system, since the drifters explore the 
circulation in the upper layer of the sea, within the first meters of water. 

On the basis of the previous considerations, we introduce our kinematic model for the 
Lagrangian dynamics. Under the incompressibility hypothesis we write a,2 — D velocity field 
in terms of a stream function: 

u = -—- and V = —-. (12) 

ay ox 

Let us write our stream function as a sum of three terms: 

^{x, y, t) = ^o(a:, y) + ^i{x, y, t) + ^2(2;, y, t) (13) 

defined as follows: 

C 
^o{x, y) = ^- [-sin{ko{y + tt)) + cos{ko{x + 2n))] (14) 

C 

\E'i(a;,|/, t) = -— ■ sin(ki{x + eisin{ujit)))sin(ki(ii + tisiniuit + 0i))) (15) 

C2 
\E'2(x, y,t) = -— ■ sin{k2{x + e2sin{u2t)))sin{k2{y + e2sin{co2t + ^2))) (16) 

k2 
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where ki = 27r/Ai, for i = 0, 1, 2, the Aj are the wavelengths of the spatial structure of the 
flow; analogously Uj = 27r /Tj, for j = 1, 2, and the Tj are the periods of the perturbations. 
In the non-dimensional expression of the equations, the units of length and time have been 
set to 200 km and 5 days, respectively. The choice of the values of the parameters is 
discussed below. 

The stationary term ^o defines the boundary large scale circulation with positive vortic- 
ity. ^1 contains the two cyclonic gyres and it is explicitly time-dependent through a periodic 
perturbation of the streamlines. The term \1'2 gives the motion over scales smaller than the 
size of the large gyres and it is time-dependent as well. A plot of the ^-isolines at fixed time 
is shown in Fig. 2. The actual basin is the inner region with negative \E' values and the zero 
isoline is taken as a dynamical barrier which defines the boundary of the domain. 

The main difference with reality is that the model domain is strictly a closed basin 
whereas the Adriatic Sea communicates with the rest of the Mediterranean through the 
Otranto Strait. That is not crucial as long as we observe the two evolutions, of experimental 
and model trajectories, within time scales smaller than the mean exit time from the sea, 
typically of the order of a few months. Furthermore, the presence of the quasi-steady cyclonic 
coastal current is compatible with the interplay between the Po river southward inflow at 
the north-western side and the Otranto channel northward inflow at the south-eastern side 
of the sea. 

The non-stationarity of the stream function is a necessary feature of a 2 — Z^ velocity field 
in order to have Lagrangian chaos and mixing properties, that is, so that a fluid particle will 
visit any portion of the domain after a sufficiently long interval of time. 

We have chosen the parameters as follows. The velocity scales Co, Ci and C2 are all 
equal to 1, which, in physical dimensions, corresponds to ~ 0.5 m/s. The wave numbers 
/cq, ki and k2 are fixed at 1/2, 1 and An, respectively. In Fig. 2a,b we can see two snapshots 
of the streamlines at fixed time. The length scales of the model Eulerian structures are of 
~ 1000 km (coastal current), ~ 200 km, (gyres) and ~ 50 km, (eddies). The typical 
recirculation times, for gyres and eddies, turn out to be of the order of 1 month and a few 
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days, respectively. 

As regards to the time-dependent terms in the stream function, the pulsations are o^i = 1 
and UJ2 = 2tt, which determine oscillations of the two large-scale vortices over a period 
Ti ~ 30 days and oscillations of the small-scale vortices over a period T2 ~ 5 days; the 
respective oscillation amphtudes are ei = 7r/5 and £2 = ci/lO which correspond to ~ 100 km 
and ~ 10 km. 

The choice of the phase factors, 0i and 02, determines how much the vortex pattern 
changes during a perturbation period. We have chosen to set both 0i and 02 to tt/4 rad. 

This choice of the parameters for the time-dependent terms in the stream function is only 
supposed to be physically reasonable, for the experimental data give us limited information 
about the time variability of the Eulerian structures. 

The chaotic advection (Ottino, 1989; Crisanti et al., 1991), occurring in our model, makes 
an ensemble of initially close trajectories spread apart from one another, until the size of 
the mean relative displacement reaches a saturation value corresponding to the finite length 
scale of the domain. 

The scale-dependent degree of chaos is given by the FSLE. Because of the relatively 
sharp separation between large and small scales in the model, we expect X{6) to display a 
step-like behavior with two plateaus, one for each characteristic time, and a cut-off at scales 
comparable with the size of the domain. In the limit of small perturbations, the FSLE gives 
an estimate of the MLE of the system. 

The LSF, i^iS), on the other hand, is expected to be proportional to the size of the 
perturbation and to A, as discussed in the introduction. Therefore the quantity 1^(6)/ 6 is 
expected to be qualitatively proportional to \{S), in the sense that the mean slopes have to 
be compatible with each other. 

In the following section we will show results of our simulations together with the outcome 
of the data analysis. 
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IV. COMPARISON BETWEEN DATA AND MODEL 

The statistical quantities relative to the drifter trajectories have been computed accord- 
ing to the following prescription. The number of selected drifters for the analysis is 37, 
distributed in 5 different deployments in the Strait of Otranto, containing, respectively, 4, 9, 
7, 7 and 10 drifters. These are the only drifter trajectories out of the whole data set which 
are long enough to study the Lagrangian motion on basin scale. To get as high statistics as 
possible, at the price of losing information on the seasonal variability, the times of all of the 
37 drifters are measured as t — to, where to is the time of deployment. Moreover, to restrict 
the analysis only to the Adriatic basin, we impose the condition that a drifter is discarded as 
soon as its latitude goes south of 39.5 N or its longitude exceeds 19.5 E. Let us consider the 
reference frame in which the axes are aligned, respectively, with the short side, orthogonal 
to the coasts, which we call the transverse direction, and the long side, along the coasts, 
which we call the longitudinal direction. 

Before the presentation of the data analysis, let us briefly discuss the problem of flnding 
characteristic Lagrangian times. A flrst obvious candidate is 

ri" = I (17) 

Of course r| is related to small scale properties. Another characteristic time, at least if 
the diffusion is standard, is the so-called integral time scale (Taylor, 1921) 

rf = -^ r C{r)dr (18) 

< f ^ > JO 

where C{t) = X^^Li < Vi{t)vi{t -|- r) > is the Lagrangian velocity correlation function and 
< f ^ > is the velocity variance. We want to stress that it is always possible (at least in 
principle) to deflne r| while to compute t^ (the integral time scale) it is necessary to be 
in a standard diffusion case (Taylor, 1921). 

The relative dispersion curves along the two natural directions of the basin, for data and 
model trajectories, are shown in Figs. 3a and 3b. The curves from the numerical simulation 
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of the model are computed observing the spreading of a cluster of 10^ initial conditions. 
When a particle reaches the boundary (^ = 0) it is eliminated. Along with observational 
and simulation data, we plot also a straight line corresponding to a standard diffusion with 
coefficient 10^ m^ / s, (Falco et al., 2000). We discuss this comparison below. Considering the 
effective diffusion properties, one should expect that the shape of Sf(t), before the saturation 
regime, can still be affected by the action of the coherent structures. Actually, neither the 
data nor the model dispersion curves display a clear power-law behavior, and are indeed 
quite irregular. The growth of the mean square radius of a cluster of drifters appears still 
strongly affected from the details of the system, and the saturation begins no later than ~ 1 
month (~ the largest characteristic Lagrangian time). This prevents any attempt at defining 
a diffusion coefficient for the effective dispersion in this system. Although the saturation 
values are very similar, we can see that, in the intermediate range, the agreement between 
observation and simulation is not good. We point out that the trouble in reproducing the 
drifter dispersion in time does not depend much on the statistics; no matter if 37 (data) or 
10"^ (model) trajectories, the problem is that the classic relative dispersion is not the most 
suitable quantity to be measured (irregular behavior even at high statistics). 

Let us now discuss the FSLE results. The curve measured from the data has been 
averaged over the total number of pairs out of 37 trajectories (~ 700), under the condition 
that the evolution of the distance between two drifters is no longer followed when any of 
the two exits the Adriatic basin (see above). In Fig. 4 the FSLE's for data and model are 
plotted. Phenomenologically, fluid particle motion is expected to be faster at small scales 
and slower at large scales. The decrease of X{6) at increasing 6 reflects the presence of 
several scales of motions (at least two) involved in the dynamics. In particular, looking at 
the values of A (5)^^ at the extreme points of the (5-range, we see that small-scale (mesoscale) 
dispersion has a characteristic time ~ 4 days and the large-scale (gyre scale) dispersion has 
a characteristic time ~ 1 month. The ratio between gyre scale and mesoscale is of the same 
order as the ratio between the inverse of their respective characteristic times (~ 10), so the 
slope of X{6) at intermediate scales is about —1. The fact that the slope is larger than —2 
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indicates that relative dispersion is faster than standard diffusion up to sub-basin scales, i.e. 
Lagrangian correlations are non-vanishing because of coherent structures. It is interesting 
to compare this Lagrangian technique of measuring the effective Lagrangian dispersion on 
finite scales to the more traditional technique of extracting a (standard) diffusivity parameter 
from the reconstruction of the small-scale anomalies in the velocity field (Falco et al., 2000). 
Estimates of the zonal and meridional diffusivity in Falco et al. (2000) are of the order of 
~ 10'^ m^/s and are compatible with the value of the effective finite-scale diffusive coefficient 
given by the FSLE, defined as X{6) ■ 6"^, computed at 5 = 20 km (~ the mesoscale). 

The FSLE computed in the numerical simulations shows two plateaus, one at small scales 
and the other at large scales, describing a system with two characteristic time scales, and 
presents the same behavior, both quahtative and quantitative, as the FSLE computed for 
the drifter trajectories. 

It is worth noting that it is much simpler for the model to reproduce, even quantitatively, 
the relation between characteristic times and scales of the drifter dynamics (FSLE) rather 
than the behavior of the relative dispersion in time. 

The LSV, in Fig. 5, shows that the behavior of i^iS), the mean velocity difference 
between two particle trajectories at varying of the scale of the separation, is compatible 
with the behavior of the FSLE as expected by dimensional arguments, i.e. v{5)/5, the LSF 
divided by the scale at which it is computed, has the same slope as A (5). This accounts for 
the robustness of the information given by the finite-scale analysis. 

We see the theoretical predictions of FSLE and LSV are fairly well comparable with the 
corresponding quantities observed from the data, if we consider the relatively simple model 
which we used for the numerical simulations. Of course, an agreement exists because of the 
appropriate choice of the parameters of the model, capable to reproduce the correct relation 
between scales of motion and characteristic times, and because large-scale (~ sub-basin 
scales) Lagrangian dispersion is weakly dependent on the small-scale (~ mesoscale) details 
of the velocity field. 
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V. DISCUSSION AND CONCLUSIONS 

In this paper we have analyzed an experimental data set recorded from Lagrangian 
surface drifters deployed in the Adriatic sea. The data span the period from December 1994 
to March 1996, during which five sets of drifters were released at different times in the vicinity 
of the same point on the eastern side of the Otranto Strait. Adopting a technique borrowed 
from the theory of dynamical systems, we studied the Lagrangian transport properties by 
measuring relative dispersions, Sf, finite-scale Lyapunov exponents A(5), and Lagrangian 
structure function i^iS). Relative dispersion as function of time does not provide much 
information, but an idea of the size of the domain where saturation sets in at long times. 
The behavior of Sf looks quite irregular and this is due not to poor statistics but rather 
to intrinsic reasons. In contrast, the results obtained with the FSLE, i.e. dispersion rates 
at different scales of motion, give a more useful description of the properties of the drifter 
spreading. In particular, X{6) detects the characteristic times associated to the Eulerian 
characteristic lengths of the system. 

We have also introduced a simple chaotic model of the Lagrangian evolution and com- 
pared it with the observations. In our point of view, the actual meaning of the chaotic model, 
in relation to the behavior of the drifters, is not that of a best-fitting model. We do not claim 
that the quite difficult task of modeling the marine surface circulation driven by wind forcing 
can be exploited by a simple dynamical system. But a simple dynamical system can give sat- 
isfactory results if we are interested in large-scale properties of Lagrangian dispersion, since 
they depend much on the topology of the velocity field and only weakly on the small-scale 
details of the velocity structures. In this respect, chaotic advection, very likely present in 
every geophysical fiuid fiow, is crucial for what concerns tracer dispersion, since it can easily 
overwhelm the effects of small-scale turbulent motions on large-scale transport (Crisanti et 
al., 1991). Even when a standard diffusivity parameter can be computed from the variance 
and the self-correlation time of the Lagrangian velocity, its relevance for reproducing the 
effective dispersion on finite scales, in presence of coherent structures, is questionable. In 
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addition, practical difficulties arising from both finite resolution and boundary effects sug- 
gest a revision of the analysis techniques to be used for studying Lagrangian motion on 
finite scales, i.e. in non-asymptotic conditions. Considering that, generally, there is more 
physical information in a scale-dependent indicator (A(5)) rather than in a time function 
{S^(t)), we come to the conclusion that the FSLE is a more appropriate tool of investigation 
of finite-scale transport properties. It is important to remark (Aurell et al, 1996; Artale 
et al., 1997; Boffetta et al., 2000) that, in realistic cases, X{6) is not just another way to 
look at S'^it) vs t, in particular it is not true that X{6) behaves like {dlnS'^(t)/dt)s'2=s'2. 
This is because X{6) is a quantity which characterizes Lagrangian properties at the scale 
(5 in a non ambiguous way. On the contrary, S'^{t) can depend strongly on 5*^(0), so that, 
in non-asymptotic conditions, it is relatively easy to get erroneous conclusions only looking 
at the shape of S^{t). In fact, at a given time, relative dispersion inside a sub-cluster of 
drifters can be rather different from other sub-clusters, e.g. because of fluctuations in the 
cross-over time between exponential and diffusive regimes. Therefore, when performing an 
average over the whole set of trajectories, one may obtain a quite spurious and inconclusive 
behavior. On the other hand, we have seen that the analysis in terms of FSLE (and LSF), 
studying the transport properties at a given spatial scale, rather than at a given time, can 
provide more reliable information on the relative dispersion of tracers. 
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FIGURE CAPTIONS 

FIGURE 1: Plot of the 37 drifter trajectories in the Adriatic Sea used for the data analysis. 
The longitude east and latitude north coordinates are in degrees. The drifters were 
deployed in the eastern side of the Otranto strait. 

FIGURE 2: Model stream function isolines at a) t = and b) t = Ti/2, with Ti ~ 30 days. 
The boundary of the domain is the zero isoline. The coordinates (x, y) are in km. 

FIGURE 3: Relative dispersion curves, for data (plus symbol) and model (dashed) trajecto- 
ries, along the two natural directions in the basin geometry: a) transverse component, 
b) longitudinal component. The time is measured in days and the dispersion in km'^. 
The experimental curve is computed over the 37 drifter trajectories; the model curve is 
computed over a cluster of 10^ particles, initially placed at the border of the southern 
gyre with a mean square displacement of ~ 50 km"^. The straight line with slope 1 
has been plotted for comparison with a standard diffusive scaling with a corresponding 
diffusion coefficient ~ 10^ m'^/s, typical of marine turbulent motions. 

FIGURE 4: Finite-Scale Lyapunov Exponent for data (continuous line) and model (dashed 
line) trajectories. 6 is in km and X{6) is in day^^. The experimental FSLE is computed 
over all the pairs of trajectories out of 37 drifters; the FSLE from the model is averaged 
over 10^ simulations. The simulated \{5) has a step- like behavior with one plateau at 
small scales (< 50 km) and one at basin scales (> 100 km.), corresponding to doubling 
times of ~ 3 days and ~ 30 days, respectively. 

FIGURE 5: Lagrangian Structure Function for data (continuous line) and model (dashed 
line) trajectories. 5 is in km and v{5)/5 is in day^^. The experimental LSF is computed 
over all the pairs of drifters; the LSF from the model is averaged over 10^ simulations. 
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